Relativistic Nucleus-Nucleus Collisions without Hydrodynamics 
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The partition function of nonequilibrium distribution which we recently obtained 
arXiv:0802.0259 in the framework of the maximum isotropization model (MIM) is exploited to 
extract physical information from experimental data on the proton rapidity and transverse mass 
distributions. We propose to partition all interacting nucleons into ensembles in accordance with 
the number of collisions. We analyze experimental rapidity distribution and get the number of 
particles in every collision ensemble. We argue that even a large number of effective nucleon 
collisions cannot lead to thermalization of nucleon system; the thermal source which describes 
the proton distribution in central rapidity region arises as a result of fast thermalization of the 
parton degrees of freedom. The obtained number of nucleons which corresponds to the thermal 
contribution is treated as a "nucleon power" of the created quark-gluon plasma in a particular 
experiment. 

PACS numbers: 25.75.-q, 25.75.Gz, 12.38.Mh 



The main goal of the investigations of the collisions 
of relativistic nuclei is extraction of a physical infor- 
mation about nuclear matter and its constituents. It 
is a matter of fact that we can get know more about 
quarks and gluons (constituents) just under extreme 
conditions, i.e. at high densities and temperatures. 
During last two decades one of the celebrated tools on 
the way of the theoretical investigations of this subject 
was relativistic hydrodynamics (RHD) which started 
to be applied to elementary particle physics from the 
famous Landau's paper 

Applying RHD one can partially describe experi- 
mental data and get know that the matter created in 
relativistic nucleus collisions (RNC) can be regarded 
on some stage of evolution as a continues one, i.e. as a 
liquid. Moreover, as was discovered in BNL, it can be 
regarded even as a perfect fluid [2| which consistent 
with a description of the created quark gluon plasma 
(QGP). The main physical quantities which can be 
extracted from experimental data exploiting the RHD 
approach are the collective (hydrodynamical) velocity 
and the elliptic flow parameter v2 of the fireball ex- 
pansion. Unfortunately, all other physical information 
is hidden in sophisticated numerical codes which solve 
Euler hydrodynamic equations of motion. 

In the present letter we propose approach to de- 
scription of relativistic heavy-ion collisions which al- 
lows to extract the physical information from exper- 
imental data on the basis of transparent analytical 
model. 

Maximum Isotropization Model. We separate 
all amount of registered nucleons into groups (ensem- 
bles) in accordance with a number of collisions, M, 
which every nucleon from a particular ensemble has 
went through. The first collision ensemble is created 
by the nucleons which take part just in one collision 
only, then M = 1, the second ensemble is created by 
the nucleons which take part in two collisions only, 
M = 2, and so on. Every ensemble contributes to 



momentum single-particle distribution function which 
can be written as 
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where the coefficients C{M) and C t h e rm reflect the 
weight of the partial contribution from every M-th 
collision ensemble and thermal distribution, respec- 
tively, to the three-dimensional momentum spectrum, 
M = corresponds to spectator particles which are 
not taking into account in this distribution. For the 
sake of simplicity we consider collision of the identi- 
cal nuclei. The expression for the partial distribution 
functions D M (p) was derived in Ref. [H, in the c.m.s. 
of colliding nuclei it reads 
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where in the Cartesian coordinate system a = 
3/2Mpj 2 nax , and in the spheric system a = 5/2Alp 2 nax , 

Z M {p) = [ d3 A . e -/3^-a[/=i + (fe,-feo,) 2 ] ^ (3) 



(3 = 1 /T is the inverse temperature, lu p = \J m 2 + p 2 , 
p^ = (jPxiPy) and Oz is the collision axis. It is un- 
derstood from ^ and (J3J) that the quantity zm(P) 
plays the role of the canonical single-particle parti- 
tion function of the M-th collision ensemble. In some 
sense particular collision ensemble M can be regarded 
as a many-particle system frozen at some stage of evo- 
lution on the way to thermal equilibrium (M cx time). 

The thermal distribution reads, -Dthcrm(p) = 
exp [-/3cj p ]/z thcrm (/3) with 2 t herm(/3) = Jd 3 ke~ /3u>k . 
Note, we separate in ([1} the thermal contribution due 
to its specific role. It would seem the contribution 
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-Dtherm(p) should appear in |T]) as a term of the ex- 
pansion with respect to partial contributions, Dm(p), 
when the number of collisions is big enough, i.e. when 
M max — > oo. Meanwhile, because of finite life time of 
the system (fireball) and hence finite number of elastic 
and inelastic collisions of nucleons this limit regime of 
hadron dynamics, M — > oo, is not achieved and M max 
is finite. However, we include in ^ a thermal source 
because it has another nature. We will return later to 
the discussion of this matter. 

There are two additional quantities in l[2" |) -([3" ]) . k^ z 
and p max , which are the external parameters deter- 
mined by the particular experimental conditions. The 
values iko z are the initial momenta of nucleons in 
c.m.s. Indeed, due to the specifics of heavy-ion colli- 
sions we know exactly the initial momenta of the nu- 
cleons in both colliding nuclei. Two Gaussians in the 
brackets on the r.h.s. of reflect a smearing around 
initial momenta which is due to collisions of nucleons 
and were obtained with the help of the saddle-point 
approximation. Note, for M = 1, 2, 3 this approxima- 
tion is not used. 

Under the notion "collision" we mean elastic rescat- 
tering as well as inelastic scattering (reactions) , which 
include a creation of secondary particles. In the trans- 
verse direction both nuclei have the same zero ini- 
tial momentum. Then, for both nuclei the smearing 
around this value is appeared in @ as the common 
factor, exp [— ap 2 j_] . The covariance of the Gaussian 
depends on the number of collisions M and the max- 
imally allowed transferred momentum. 

The rapidity distribution was obtained after inte- 
gration of {T]) with respect to the nucleon transverse 
mass, m± = (m 2 N + p 2 ^) 1 ^ 2 , where itin is the nu- 
cleon mass and rapidity, y, is defined as tanhy = 
P z /lo p . With respect to new variables one gets, d 3 p = 
d(j>u)p m± dm± dy. 

As a first step of our approach we consider a central 
collision of identical nuclei and we assume an azimuth 
symmetry of radiation of the particles. Rapidity spec- 
trum of registered particles looks like 
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C{M) IPM(V) + Cthorm ^therm(2/) , (4) 



M=M„ 



where 



ip M (y) = 27r cosh y / dm^ m\ D M {m± , y) . (5) 

To obtain ^>therm(2/) we put -Dtherm in place of Dm on 
the r.h.s. of |5]). Double differential spectrum which 
depends on the transverse mass is obtained from ([TJ) 
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Usually the mode M = 1 does not give contribution to 
the particular experimental rapidity window. In this 
case we can set C(l) ~ and start summation in (j4| 
from M m ; n = 2. 

Extraction of the physical information from 
experimental data. With making use of the ra- 
pidity distribution ^ we fit the experimental data 
on the rapidity distribution of net protons which were 
measured at the CERN SPS (NA49 Collaboration) 
The slope parameter (3 was first extracted from 
double differential yield for protons with the use of 
the thermal distribution. The proton data is remark- 
able in that sense that we know exactly the initial 
momentum, kg z , of every nucleon. The fit was car- 
ried out with a help of the program MINUIT, vari- 
able parameters are coefficients C(M) and C t h e rm 
and parameter, which confines the momentum space, 
Pmax- The values of the obtained parameters for 
T = 1/13 = 248 MeV are shown in Table 1. All evalu- 
ations are carried out in the c.m.s. of colliding nuclei 
with use of the spheric coordinate system. 

Table 1 
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1.275 
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The obtained theoretical curves together with exper- 
imental data are depicted in Fig. [1] Broken curves 
(see upper panel) marked by the numbers M and 
solid thick curve (blue in on-line presentation) rep- 
resent the partial contributions from every ensemble, 
£7(M) -(pM(y)j an d complete theoretical proton rapid- 
ity distribution, respectively. The thermal contribu- 
tion is represented by central bell-like dashed curve 
(red in on-line presentation). 

The integral on the r.h.s. of Eq. ([3]) which gives 
rise to single-particle partition function zm((3) is de- 
fined in the rapidity range [— l^ m ,i^m]; where F cm = 
^bcam/2. Then, the functions ifM(y) are normalized 
to unity in the same range. If one integrates Eq. (0} 
with respect to rapidity in this range it is easy to find 
that a result of integration on the r.h.s. equals to 
the sum of all coefficients C(M) plus Ctherm- At the 
same time the value of this integral equals the area 
under the "rapidity" curve (solid, thick blue curve) 
in FigfT] or to the total number of participated pro- 
tons which would be registered in case if the total ra- 
pidity window [-Fcni! ^cm] is allowed experimentally: 
^tot) = ^ c ^ + Cthcim ^ For instancej for N A49 

experimental data we obtain Np tot ^ — 151. 

So, every coefficient £7(M) tells us how many pro- 
tons undergo M effective collisions or what is the pop- 
ularity of every collision ensemble. For instance, en- 
semble of the protons which participated in nine ef- 
fective collisions, M = 9, consists of 10 protons, i.e. 
£7(9) ~ 10. What is very important, we learn from 
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SPS(NA49): Pb+Pb, E kin =158AGeV 
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FIG. 1: Upper panel: The result of the fit (thick solid 
curve, blue in on-line presentation) to experimental data 
U on rapidity distribution (see Eq. ©). Broken curves 
marked by the numbers M represent the partial contri- 
butions from every collision ensemble, C(M) • tfM{y), the 
thermal contribution is represented by central Gaussian- 
like dashed curve (red in on-line presentation). 
Lower panel: The thick solid curve (blue in on-line pre- 
sentation) represents the evaluation of the mj_-spectrum 
obtained in accordance with Eq. (J5J) were we use the same 
values of the coefficients C{M) which were obtained as 
a result of the fit of the dN/dy data. Experimental data 
points are from [j. Broken curves marked by the numbers 
M represent the partial contributions from every collision 
ensemble, 2nm± coshj/ • C(M)DM{m±,y). 



this expansion that Ctherm — 18, it means that ap- 
proximately eighteen protons come from a thermal 
source what makes up 12% (twelve percent) of all par- 
ticipated protons. 

Now we would like to draw attention to the ensem- 
ble with maximal collisions, M = 12. This value was 
determined from UrQMD ;5| evaluation of the mean 
maximum number of the effective nucleon collisions. 
It turns out, as it is seen from FigfT] (upper panel), 



that the partial function <fi2(y) does not "fill in" suc- 
cessfully the central rapidity region. That is why the 
presence of the thermal function Ctherm (y) (Gaussian 
like curve in the center of FigfT]) is necessary in the 
expansion ([4]). Even this big number of effective colli- 
sions, M max = 12, cannot give rise to a source which 
is compared in the central rapidity region with a ther- 
mal one. It is the main reason why the thermal source 
should be presented in the expansion ((4]). 

Note, in the case of finite and small number of 
experimental points the set of functions, tfM(y)i is 
overcomplete. To choose a unique configuration of 
the variable parameters we use the maximum entropy 
method @- 

In this analysis we are coming to one of the main 
conclusions, which can be derived from our model: 
The thermal source has absolutely different nature 
of origination, it cannot be created just due to the 
hadron reactions of nucleons which result in random- 
ization and subsequent isotropization of the nucleon 
momentum. The thermal source can emerge as a re- 
sult of appearance "at once" of many new degrees of 
freedom. We know just one candidate to this role, it 
is the quark-gluon plasma, for instance, its creation 
can occur in collision of nucleons in the presence of 
a dense medium, N + N — > n g + n q . Then, a many- 
parton system, which emerges in the collision, consists 
of n g ^> 1 gluons and n q ^> 1 quarks. All momenta 
of quarks and gluons can be regarded from the very 
beginning as random ones and thermalization of the 
system occur during a time span T t hcrm = 0.6 fm/c 
[8| . Hence, the protons which come from the thermal 
source indicate the presence of the QGP in the fire- 
ball and we can determine a power of the QGP by the 
number of protons outcoming from the thermal source 
or by the value of Ctherm- 

Actually, the total number of nucleons which ap- 
pear as a result of hadronization of the QGP can be 
then evaluated with accounting for isotope composi- 



tion of the colliding nuclei: N 



(QGP) 
N 



c, 



thcrr 



.4 
Z 



For 



instance, in the experiment under consideration we 

find Ctherm — 18, and then iV^ GP ^ ~ 46, i.e. ap- 
proximately 46 nucleons were created by the QGP or 
by several QGP drops. This makes up 12% from a 
total number of net nucleons which are the partic- 



ipants of the collision, N 



(participants) 
N 



382. Then, 



we estimate a "nucleon power" of the QGP, P qgp = 
A r ^ GP ' ) /A r ^ artlclpants \ which was created in nucleus- 
nucleus collision 



P. 



12% 



o m 



For instance, it turns out 
Pb+Pb collisions (SPS) at E kin 



that 



158A GeV. 

The same analysis was carried out for proton dis- 
tribution from 11. 6 A GeV/c Au + Au collisions mea- 
sured by the E802 Collaboration 0. The fit to the 
experimental data (0-3% centrality) allows to extract 
the values of parameters which are shown in Table 
2. For this experiment UrQMD 0] evaluation gives 
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the mean maximum number of the effective nucleon 
collisions M max = f3. The width of the rapidity 
window in this experiment avoids the necessity to 
take into account the collision ensemble M — 1 too. 

Table 2 



AGS(E802): Au+Au, p, =11.6A GeV/c 
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C(3) 


(7(4) 


(7(5) 
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(7(9) 


12.2 


14.5 


6.5 


5.1 


5.1 


5.8 
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8.2 



(7(10) 


C(ll) 


£7(12) 


C(13) 


Cthcrm 


Pmax (GeV/c) 


9.7 


11.4 


13.0 


14.7 


37.9 


0.724 



Theoretical curves and experimental data are depicted 
in Fig. [2J Notations and marks have the same mean- 
ing as in the previous figure. In the lower panel 
the solid curves represent the evaluation of the mj_- 
spectra in different rapidity windows. Remind, these 
curves are obtained without additional fitting of the 
data. We just use the values of parameters from the 
Table 2. Meanwhile, on the first step, before fitting 
the rapidity distribution, we estimate with the help 
of the thermal distribution the slope parameter (3 and 
find T = 280 MeV. 

We can estimate as well the nucleon power of the 
produced QGP in experiment Au + Au at 11. 6 A 
GeV/c (0-3% centrality). In this Case Ctherm — 38 



and N, 



(tot) 



151, hence we obtain P c 



qgp 



25%. 



Summary and discussion. In the proposed 
Maximum Isotropization Model [|[ the maximum 
number of collisions (reactions), Af max , assumed to 
be finite and determined by the nuclear number A, 
initial energy and centrality. With the help of the 
UrQMD transport model [5[ it was found that for 
SPS (Pb+Pb, 158A GeV) conditions Q, M max = 12, 
and for AGS (Au+Au, 11.6A GeV/c, 0-3% centrality) 
conditions [6j, M max = 13. Utilizing thermal distribu- 
tion we extract a slope parameter from experimental 
data on the proton mj_-spcctra: for SPS conditions 
0, T = 248 MeV, and for AGS conditions (0- 
3% centrality), T = 280 MeV. It is evidently seen 
from Fig. [T] (lower panel) that the mj_-spectrum is 
mainly determined by the thermal component, and in 
any case the slope of the partial contribution, marked 
by M, is approximately the same as of thermal dis- 
tribution. Exactly of that reason the m_L-spectrum is 
low informative about collision ensembles or the infor- 
mation about rescattered nucleons almost lost in this 
presentation. On the other hand, the m_L-spectrum 
as a trigger gives possibility to extract a value of the 
slope parameter. 

Next, we made the fit of experimental data [3,0 on 
the rapidity distribution of the net protons and ob- 
tained the set of coefficients C(M) (see Tables 1, 2) 
which are nothing more as an absolute number of 
protons in every collision ensemble. Note, the pro- 
ton data is interesting first of all because we know 
an exact value of the initial nucleon momentum. As 
a matter of fact, the partial expansion, dN/dy = 
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FIG. 2: Upper panel: The result of the fit (thick solid 
curve, blue in on-line presentation) to experimental data 
[6( on rapidity distribution (see Eq. Q). Broken curves 
marked by the numbers M represent the partial contri- 
butions from every collision ensemble, C(M) ■ tp M (y), the 
thermal contribution is represented by central Gaussian- 
like dashed curve (red in on-line presentation). 
Lower panel: The solid curves represent the evaluation of 
the m±-spectra obtained in accordance with Eq. ([6]) were 
we use the same values of the coefficients C(M) which 
were obtained as a result of the fit of the dN/dy data. 
Experimental data points are from [6|. 



tained coefficients C(M) from Tables 1, 2, makes a 
good description of the experimental data on rapid- 
ity distribution, except the central rapidity region. 
It means that the central rapidity region cannot be 
described just by finite number of nucleon rescatter- 
ings (hadron reactions). Then, we are forced to take 
into account also the thermal source, which evidently 
has a different nature. We assume that this source is 
a thermalized multi-parton system (QGP) @ which 
through hadronization process emits totally thermal- 
ized nucleons. The knowledge of the number of pro- 
tons, Cthorm, which come from the QGP, gives us 
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a possibility to evaluate the "nucleoli power" of the 
QGP, fqgp, created in the particular experiment on 
nucleus-nucleus collision. We find that for SPS con- 
ditions [H, -Pqgp — 12%, and for AGS conditions [(J 
(0-3% centrality), P qgp 25%. So, following the pro- 
posed criterium we can claim that QGP (as a nucleoli 
source) was created not only at SPS energies but 
it was also created, even more powerful with respect 
to nucleons, in the central collisions at AGS energies. 
Meanwhile, in accordance with UrQMD estimations 
the number of pions created in hadron reactions at 
the SPS [H and AGS @ energies are approximately 
the same. Hence, the number of pions created by the 
thermal source at the SPS is much bigger than the 
number of pions created by the thermal source at the 



AGS. From that we can conclude that "pion power" of 
QGP created in nucleus-nucleus collisions at the SPS 
up to one order higher than that one created at the 
AGS. 

All this leaves us with the continued challenge of 
applying the model to other experiments and prob- 
lems. 
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